Polarizability and dynamic structure factor of the one-dimensional Bose gas 
near the Tonks- Girardeau limit at finite temperatures 
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Correlation functions related to the dynamic density response of the one-dimensional Bose gas 
in the model of Lieb and Liniger are calculated. An exact Bose-Fermi mapping is used to work in 
a fermionic representation with a pseudopotential Hamiltonian. The Hartree-Fock and generalized 
random phase approximations are derived and the dynamic polarizability is calculated. The results 
are valid to first order in I/7 where 7 is Lieb-Liniger coupling parameter. Approximations for 
the dynamic and static structure factor at finite temperature are presented. The results preclude 
superfluidity at any finite temperature in the large-7 regime due to the Landau criterion. Due to 
the exact Bose-Fermi duality, the results apply for spinless fermions with weak p-wave interactions 
as well as for strongly interacting bosons. 
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I. INTRODUCTION 

The one-dimensional (ID) Bose gas with point inter- 
actions is an abstract and simple, yet nontrivial model of 
an interacting many-body system. The recent progress 
in the trapping, cooling, and manipulation of atoms have 
made it possible to test theoretical predictions experi- 
mentally and have thus led to a revival of interest in this 
model. 

Exact solutions were described by Lieb and Lin- 
iger who found that the physics of the homogeneous 
Bose gas is governed by the single dimensionless parame- 
ter 7 = g^ra/ \h 2 n). Here n is the linear particle density, 
m is the mass, and c/b is the strength of the short-range 
interaction between particles For small 7 we have a 
gas of weakly interacting bosons. Although there is no 
Bose-condensation in one dimension even at zero tem- 
perature Q , many properties of the gas are reminiscent 
of Bose-Einstein condensates with Bogoliubov perturba- 
tion theory being valid and even superfluid properties 
were predicted [HIE For large 7, however, the sys- 
tem crosses over into a strongly interacting regime and 
at infinite 7 we obtain the Tonks-Girardeau (TG) gas of 
impenetrable bosons In this regime, the strongly re- 
pulsive short-range interaction has the same effect as the 
Pauli principle for fermions. Many properties like the ex- 
citation spectrum become that of a free Fermi gas and, 
indeed, the model maps one-to-one to a gas of noninter- 
acting spinless fermions. 

Experiments have recently probed the crossover to the 
strongly-correlated TG regime by increasing interactions 
up to values of 7 w 5.5 and to effective values of 
7 c ff sa 200 in an optical lattice |K| ■ The momentum dis- 
tribution [Tol |. density profiles Rj, and low-energy com- 
pressional excitation modes have been the focus of 
the experimental studies. In a recent experiment |l2j |. 
the zero momentum excitations of a ID Bose gas in an 
optical lattice have been measured by Bragg scattering, 



a technique that could also be used to measure the dy- 
namic structure factor (DSF) S(q,iu), which is calculated 
in this paper [T^| . 

The theoretical description of the Lieb-Liniger model 
is not complete. Although the exact wavefunctions, 
the excitation spectrum, and the thermodynamic prop- 
erties are known for arbitrary values of the cou- 
pling constant 7, it is notoriously difficult to calcu- 
late the correlation functions. Many results in limit- 
ing cases are summarized in the book |l5j . but the full 
problem is not yet solved. Recently, progress has been 
made on the large-distance and long-time asymptotics of 
single-particle correlation function s 11^11^1 a nd on time- 
independent correlation functions [SHliS EJ HI HI 
l24j . In the general case and for time-dependent corre- 
lation functions, a wealth of information is available for 
small 7 where Bogoliubov perturbation theory can be ap- 
plied as well as for the TG gas at 7 = 00. However, the 
strongly-interacting regime with large but finite 7 was 
hardly accessible as a systematic expansion in 7 _1 was 
lacking. 

In ID there is a duality between interacting Bose and 
Fermi man y-b ody systems. A couple of recent works 
[H HI ElSl Hi El already pointed out that the ex- 
act Bose-Fermi mapping that Girardeau used to solve 
the case of 7 = 00 Q can be extended to the case of 
finite interaction 7. Thus, a model system of interacting 
fermions can be constructed for which the energy spec- 
trum and associated wavefunctions are in a one-to-one 
correspondence with the Lieb-Liniger solutions. Our ap- 
proach makes use of the same Bose-Fermi duality, how- 
ever, the motivation is to derive new results for the ID 
Bose gas. We can calculate correlation functions of the 
strongly-interacting Bose gas by solving the equivalent 
interacting Fermi problem in the regime where its inter- 
actions are small. In our previous paper [3lj we calcu- 
lated the DSF for the Lieb-Liniger model for large 7 at 
zero temperature and related it to the Landau criterion of 
superfluidity. The DSF S(q,u>) holds information about 
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the strength or the excitability of excitations with mo- 
mentum hq and energy fko and thus may indicate decay 
routes of possible supcrcurrents. Although the TG gas 
is not supcrfluid due to low-energy umklapp excitations 
near 10 = and q = 2wn, a crossover to superfluid be- 
haviour for finite 7 is possible if the umklapp excitations 
are suppressed. This possibility is not precluded by the 
results of Rcf. jp. 

The purpose of this paper is to calculate so far un- 
known correlation properties of the ID Bose gas at finite 
temperatures in the strongly-interacting regime. We cal- 
culate the dynamic density-density response, the DSF, 
and the static structure factor of the Lieb-Liniger gas 
with the fermionic random-phase approximation (RPA), 
extending our previous results |3 1| to the case of fi- 
nite temperatures. Although our calculations are non- 
perturbative, we obtain the first order term in the expan- 
sion in 7 _1 for comparison. In particular we find that the 
ID Bose gas at finite temperature cannot be superfluid 
due to the finite probability of umklapp excitations. We 
also present an extended discussion of the validity of the 
pscudopotcntial approach and give a critical scrutiny of 
the limits of applicability of our approach. 

The structure of this paper is as follows. In the next 
section we consider the exact Bose-Fermi mapping for 
finite values of the interaction strength and discuss the 
use of a fermionic pseudopotential. In Sec. II I II we de- 
rive a Hartree-Fock (HF) approximation for the fermions. 
The generalized Random-Phase Approximation (RPA) 
is derived in Sec. IIVI as a linearized time-dependent HF 
scheme. Analytic expressions for the polarizability, the 
DSF, and the static structure factor arc analyzed and 
discussed. 



II. FERMIONIC PSEUDOPOTENTIAL 

We consider the system of N interacting bosons of mass 
m in ID described by the Hamiltonian 



Hr- 



N 

E 



H 2 d 2 
2m dx 2 



g-B*^2S(xi - xj). 

i<j 



(1) 

This model extends the Lieb-Liniger model 0, 0] to in- 
clude an external potential V ex t(x)- In contrast to the 
Bethe ansatz solutions of Refs. 0, Q, our approach de- 
veloped below allows us to treat the effects of external 
potentials, which are important in the context of experi- 
mental realizations. 

We will now map the model Q onto an equivalent 
Fermi system and discuss the appropriate pseudopotcn- 
tials. As described in detail in Ref. pj, the 5-function 
point interaction in the Hamiltonian Q can be repre- 
sented as a boundary condition for the wavefunction in 
coordinate space at the points where two particles meet 
at the same position. For simplicity we will first dis- 
cuss the case of two particles and introduce center-of- 



mass R = (x\ + X2)/2 and relative x = X2 — x\ coordi- 
nates. The bosonic wavefunction has the (even) symme- 
try ip B (x,R) — ip B (—x, R). Due to this symmetry, the 
effect of the 8- interaction on ip B (x, R) can be formulated 
as a single boundary condition for x — ► +0: 



lim d x ip(x,R) 
v— >+o 



lim 



o~2fi2~ 



(2) 



with i/j(x,R) = ^(xjR)- The exact Bosc-Fermi map- 
ping now takes advantage of this boundary condition 
being formulated for x > where ip B (x,R) solves 
the Schrodingcr equation. A fermionic model is de- 
fined by the same boundary condition @ and the same 
Schrodinger equation for 1 / but requiring fermionic 
(odd) symmetry. We thus find fermionic solutions ifj F 
with 



ip F (x,R) 



ip B (x,R), x > 0, 
-ip B (x,R), x < 0. 



(3) 



The fermionic symmetry together with the boundary 
conditions J2J) requires a discontinuity in the wave- 
function ip F (x : R) at x = inducing a jump of 
4fi 2 / (gBm)d x tp F (x = 0,R) in the wavefunction and a 
continuous first derivative, whereas the bosonic wave- 
function ip B (x,R) is continuous but has a discontinuous 
first derivative. In the simple limiting case of <?b —> °o 
we obtain the TG gas and i/j f is continuous. 

This generalized Bose-Fermi mapping was introduced 
by Cheon and Shigehara j2f|, who also discussed the 
straightforward generalization to the iV-particle problem 
of Eq. Q . The mapping as described above is exact and 
one-to-one for particles confined by an external potential 
V CX f If periodic boundary conditions are imposed upon 
the Bose system, they translate in the Fermi-system into 
periodic boundary for odd TV and antipcriodic boundary 
for even N [sec Eq. (37) of Ref. [2f|]. Antiperiodic bound- 
aries mean that the fermionic wavefunction changes by 
a factor of —1 whenever a particle is translated by the 
length of the periodic box L. The differences between pe- 
riodic and antiperiodic boundaries are, however, minute 
for large systems and vanish in the thermodynamic limit. 
For this reason we only consider explicitly the case of pe- 
riodic boundary conditions for the fermionic wavefunc- 
tions below. 

As a result of the Bose-Fermi mapping, the energy 
spectrum of the Bose and corresponding Fermi system 
are identical. Furthermore, all observables that are func- 
tions of the local density operators are identical in both 
systems because they involve absolute values of the wave- 
functions only and sign changes as in Eq. (|SJ) do not mat- 
ter. In particular, this includes the dynamical density- 
density correlation functions and derived quantities like 
the dynamic and static structure factors. By contrast, 
the off-diagonal parts of the one-body density matrix and 
consequently the momentum distribution show distinct 
differences in both systems |32l |. 

For our purposes it is desirable to represent the interac- 
tion in the fermionic model as an operator. In fact, it has 
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been shown by Seba that the discontinuity-introducing 
boundary condition (2} for fcrmionic symmetry defines 
a self-adjoint operator on Hilbert space Seba also 

gave an explicit construction by the zero-range limit of a 
renormalized separable operator of finite range. As a re- 
sult, one can represent the fermionic interaction in terms 
of an integral kernel pF * 



V F {x 1 ,x 2 \x' 2l x' l ) 



6'(x 1 -x 2 )5'(x' 1 -x' 2 ), 



(4) 



where the coupling constant in the fermionic representa- 
tion is defined as 



g F = 2H 4 /(m 2 g B ) 



(5) 



Due to the gap in the fcrmionic wave functions 10. we 
should be very careful defining the matrix element of 
the pseudopotcntial for two arbitrary fermionic two- 
particle wavefunctions: 



(V F |VF^ F )=-2.g F / dR[d x il> F (x,R)]*d x <p F (x,R) 



(6) 

where due to the fermionic symmetry the right and left 
limits of the first derivatives \im x ^± d x ip F (x, R) co- 
incide, even if the wavefunctions are discontinuous at 
x = 0. Representations similar to Eqs. 0} and |jf)} have 
also been given in Refs. [28|,|2?],|3(|. The fermionic Hamil- 
tonian H F takes the form of Eq. but with the inter- 
action term J2i<j V~f{ x i 

) instead of the bosonic 
(^-function interactions. If the bosonic wavefunction ip^ 
is an eigenfunction of the Hamiltonian Hb of Eq. (|TJ 
with eigenvalue , then the fermionic wavefunction ifj F 
of Eq. ® is also eigenfunction of the fermionic Hamilto- 
nian with interaction V F of Eq. JIJ corresponding to the 
same eigenvalue E F = E*. This can be verified easily 
for two particles by substituting ip F into the Schrodingcr 
equation. We thus conclude that the fcrmionic represen- 
tation with interactions 0} is exact for all values of the 
coupling constant g F . 

In addition to the formal considerations above we now 
give another, somewhat heuristic justification for the 
pseudopotential (QJ following Sen's argument [34| based 
on the Hellmann-Feynman theorem. Let us suppose that 
we know the exact eigenvalue E® of Hb together with the 
corresponding bosonic wave function if^ . Then it follows 
from the Hellmann-Feynman theorem that 



dgB 



dH F 



= / dR 



dgB 

9 2~ 
m <?B 



<P 



= J dR\?p B (x = ±Q,R)\ 2 
\d x ij F (x = ±0,R)\ 2 , (7) 



where Eqs. J3} and © were used to derive the last equal- 
ity. With the help of Eqs. JSJ) and JBJ we find dE®/dg F = 



-2(i; F \dV F /dg F W) and finally 8E*/dg F = dE F /dg F 
by the Hellmann-Feynman theorem for the eigenvalue E F 
of H F . Taking into account that Hb and H F have the 
same eigenvalue spectrum in the TG limit of g F — > 0, 
we conclude that both Hamiltonians have the same spec- 
trum also for arbitrary values of g F . 

A useful approximate representation as a pairwise pseu- 
dopotential was suggested by Sen 



V Se n{xi,X 2 ) = -g F S"(xi - X 2 ), 



(8) 



where S"(x) denotes the second derivative of the delta 
function. In contrast to the integral kernel 0}, Sen's 
pseudopotential takes the form of a local operator, which 
yields a simplification when performing analytical cal- 
culations. This pseudopotential, however, is applicable 
only for variational calculations in a variational space of 
continuous fermionic functions that vanish whenever two 
particle coordinates coincide. This is the case for Slater 
determinants that may be used to derive the Hartree- 
Fock (HF) and Random-Phase approximations (RPA) 
but not for the exact fermionic wave functions like l|2}. 
One can justify Sen's pseudopotential JH} up to first order 
+±o' in 7 _1 in the same manner as in the previous paragraph. 



For higher orders, Vs en is not correct, because its matrix 
elements contain not only the correct term, as in the r.h.s. 
of Eq. ©, but an additional term disappearing only at 
7" 1 = 0. 

The fcrmionic Hamiltonian can now be rewritten in 
terms of Fermi field operators ^>(x) and ^(x) 

h f = [ dx 9_^imm + / dxVcxt (x)^(xMx) 



dxida^dx'jdxj Vp(x\, x 2 ;x' 2 ,x' 1 ) 



x ¥(x 1 )&{x 2 )i>{x' 2 )i>(x' 1 ), 



(9) 



where V F is given by Eq. lf3|. Alternatively, the approx- 
imate pseudopotential Vsen can be employed [3f|. 

In the remainder of this paper we will study the 
fermionic model © in the HF approximation and the 
RPA. 



III. THE HARTREE-FOCK OPERATOR 

The HF approximation for the fermionic system (|5} 
is derived in the standard way by variation over Slater 
determinants. We thus expect Sen's pseudopotential (jHJ 
to be valid. Indeed, we find that the interactions (|S} and 
Q yield identical results on the HF level. 

When working at finite temperatures, it is convenient 
to introduce the HF operator as the single-particle oper- 
ator 



H = dxdx' F(x,x')¥(x)$(x') 



(10) 
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that minimizes the Gibbs-Bogoliubov inequality [36j with 
respect to F(x, x') 

n^Q + (H F - H ) . 

Here f2 = —a In Z is the grand thermodynamic potential 
with the partition function Z = Tr cxp[— j3(H F — fiN)] 
corresponding to the Hamiltonian H F . The inverse tem- 
perature P = 1 /T is introduced here and we use the units 
ks = 1 in this paper. Accordingly, VLq is associated with 
Ho, and (...) = ^ Trf • • exp[-/3(# - fxN)}}. The 
variational procedure yields 



h 2 d 2 



F(y,z) = S(y-z) 

+ / dxdx 1 [V{y,x;x',z)-V{y,x;z,x')]p( 1 \x,x l ), 



(11) 



with the one-body density matrix p^(x',x) = 
(x')^>(x))o, which should be determined in a self- 
consistent manner. The simplest way to do this is to 
work in the diagonal representation of the HF kernel 
F(x,x') = J2j £ j ( Pj( x ') l Pj( x )- The single-particle func- 
tions <fij(x) are called Hartree-Fock orbitals. This rep- 
resentation allows us to rewrite the HF Hamiltonian 
dll) [I in terms of the creation and destruction operators 
a! = J dx& (x)ifj(x) and a,j = J dx &(x)<p*(x), respec- 
tively. It takes the form Hq — Y^j £jojoj, which leads to 



where 



(aj Oj) = rijSij, 



exp[/3(£j - /ti)] + 1' 



(12) 



(13) 



is the Fermi distribution of occupation numbers, and Sij 
is the Kroneker symbol. At zero temperature, rij defines 
the Fermi step function. By using the representation 
= Ylj ^jfji 3 ^) an d Eq. (|12|l . we derive 



p (1 \x lX ') =^2rij(p* j (x)ip j (x'). 



(14) 



By substituting Eq. i|14|) and either one of the the pscu- 
dopotentials © or (0} into Eq. (|llfl . we come to the same 
local form of the HF kernel 



F 



d 2 d 
n{x)—+2V{x)i—-T{x) 



9f 



dx 2 



dx 



(15) 



which is defined as Fip(x) = J dx' F(x, x')f(x'). The 
first two terms on the right hand side come from the 



single-particle part of the Hamiltonian ©. The mean- 
field parts in the square bracket involve the local den- 
sity n(x) — p lyl \x,x) = nj<P* j(x)(fij(x) and the 

derivative densities V(x) = —iJ2j n j ( P*'j( x )'Pj( x ) an d 
T(x) = J2 3 nj[^*{x)^{x) + 2p*' j (x)ip' j (x)}. Here we de- 
fine ip' = dip/dx. The quantities V(x) and T{x) are 
reminiscent of momentum and energy densities, respec- 
tively. We find a purely local Fock operator, contrary to 
the case of Coulomb interactions where the Fock oper- 
ator F is nonlocal with a local Hartree and a nonlocal 
exchange term. 

For the homogeneous gas {V ey ± = 0), the quantum 
number j can be associated with the particle HF orbitals 
are plane waves <p q (x) — cxp(ixq)/^/L with energy and 
effective mass 



h 2 q 2 
2m* 



- 9vn{q 
m 



1 - 2g F mn/H 2 1 - 47- 



(16) 
(17) 



respectively, where we have introduce the average square 
momentum 



(18) 



over the Fermi distribution n q of Eq. i|13|) . Here, the 
sum over q runs over values q — 2irl/L, I = 0, ±1, ±2, • ■ • 
in accordance with periodic boundary conditions. The 
chemical potential and the density cannot be indepen- 
dent quantities; they arc related through 



(19) 



In the thermodynamic limit n = N/L = const, L — > 
00 all the sums over momentum become integrals: 
(1/L) J2 q ~ ¥ (2?r) / dq ■ ■ ■ . In the canonical ensemble 
only two thermodynamic parameters are independent, 
the density and temperature. Thus, we can use 7 and 
f3 as input parameters and determine the chemical po- 
tential p in a self-consistent manner from Eqs. (|16|) - (|19fl . 

At zero temperature the HF scheme admits the ana- 
lytical solutions (q 2 ) = fc F /3, p = Ek F = s F [l — 16/(37)], 
where k F = im and energy e F = fi 2 fc F /(2m) are Fermi 
wavenumber and energy of the TG gas, respectively. The 
HF approximation for the ground-state energy yields 

?} 2 7r 2 n 2 

E KF =(H F ) = N— (I-87- 1 ). (20) 

om 

It coincides with the first two terms of the largc-7 expan- 
sion of the exact ground-state energy in the Licb-Liniger 
model Q. Note that Eu F (-Ho)o = J2 P £ p-> as it should 
be in the HF scheme (see e.g. Ref. |13)- 

We now briefly discuss the stability of the HF solu- 
tion. Stability of the HF solution implies positivity of the 
isothermal compressibility of the medium (dn/ &p)t /n. 



The latter relates directly to the isothermal speed of 
sound 



VT = —(^) 

V m \dnJ t' 



(21) 



shown in Fig. ^ for various temperatures. The HF result 
at T = 



Vt = \J 1 — 87 1 

m 

= — [l-^ + Ofr" 2 )] 



(22) 



yields the correct first order expansion of Lieb's exact 
result 0. We see in Fig. ^ that the stability condition 
v\ > is broken below some critical value of the coupling 
constant 7, depending on temperature. From Eqs. Ilti|) - 
(|19|1 one can show that the critical value lies between 
7 = 8 at zero temperature and 7 = 4^/ (1 — v3) ~ 9.464 
at large temperatures. Thus, the developed HF scheme 
(and, hence, the RPA discussed below) is applicable only 
for values of the Lieb-Linigcr coupling constant of the 
order 7 > 10. 
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IV. THE RANDOM PHASE APPROXIMATION 

A. Response function 

The HF approximation permits us to calculate the lin- 
ear response of time-dependent HF. Approximations of 
the linear response functions on this level are known 
as RPA with exchange or generalized RPA [3^, |3^ ■ In 
this section we will calculate the density-density response 
function x(<Z, z) also known as dynamic polarizability. It 
is intimately related to the DSF and the time-dependent 
density-density correlation function [3^, |3^, . In order 
to define x(g, z) we consider the linear response of the 
density 

n(x, t)-n= ($t(s, 0{x, t))-n= - ^ e iqx Sn(q, t) 

1 

to an infinitesimal time-dependent external potential 

q J 71 

Here we choose e — > +0 to provide the boundary condi- 
tion T4 x t(x, t) — * when t — * —00. For q 7^ we have 

Sn(q,t)= ^2{a{_ q/2 (t)a k+q/2 (t)) 

k 

= J fi^e- iut e et 8n(q,u). (23) 

The dynamic polarizability is now defined by 

X(q, uj + ie) = -Sn(q, uj)/5V ext (q, ui) (24) 



i , 1 , 1 1 1 , 1 

0.05 0.1 0.15 0.2 

l/y 

FIG. 1: (Color online) (a) The isothermal speed of sound vt 
of Eq. 12111 versus the inverse coupling constant 7 -1 in the 
HF approximation for different temperatures. The quantity 
Do = him/m is the speed of sound of the TG gas at zero tem- 
perature. The speed of sound becomes zero at some critical 
value of 7, below which the HF solution becomes unstable (see 
discussion in Sec. IIIII 1 (b) The solid (black) line shows the 
speed of sound in the HF approximation at T — given by 
Eq. 12211 . and the dashed (violet) line shows the exact speed 
of sound in the Lieb-Liniger model for comparison. 

and obviously determines the linear density response to 
an external field. 

The polarizability can be obtained directly from the 
linearized equation of motion of the density operator in 
the time-dependent HF approximation in the standard 
way as summarized below. 

(i) With the help of the HF Hamiltonian i|10|) we write 
the equation of motion ihdp^/dt = [p^'jHo] for the 
operator = & (y , (z , t) and take its average. We 
thus derive 

ift— — 7^ Z ' ~ = J dr [F(z,x)p {1) (y, x,t) 

-F(x,y)pW(x,z,t)] (25) 

with the HF kernel F of Eq. (TTJ. 

(ii) We substitute p^(y,z,t) = po(y — z) + 
5p^\y,z,t) into Eq. (|25fl and linearize it with respect 
to 5p^ and <5Vj, x t. Here we introduce the equilibrium 
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value of the one-body density matrix p^ (y — z) = 
(1/L) n k exp[ik(y — z)] in the HF approximation with 
the HF occupation numbers n k of Eq. (|T3|) . 

(iii) In the Fourier representation of momentum and 
frequency, the obtained linearized equation becomes al- 
gebraic and takes the form 

6pW (k, q, U ) (\Y. n P [V(* - P - 9/2) -V(k-p + q/2)\ 



h 2 kq 



fku — is 



X {n k +q/2 - Tlk-q/2)- 



(26) 



Here, V{q) = g F q 2 stands for the Fourier transform of 
the potential © and 5p^ is defined by the relation 
(4-«/a(*) 6 *+«/3(*)> = J dLoe-^e^SpW^q^) 

for g / 0. We are interested in the density response 
Sn(q,u>), which is directly connected to Sp^ by 



Sn(q,uj) = ^(5p (1) (fc,g,w). 



(27) 



Because V(q) is a polynomial in q, wc can obtain 
an analytical expression for the polarizability l|24(l from 
Eqs. (|26[l and (|27|l . After somewhat lengthy but straight- 
forward calculations wc find 



x(q,z) = 

with z = lu 



(l-4 1 -^[B + D(q 7 z) X ^(q,z)} 
is and s^ + 0. Here we denote 



(28) 



B= l-4(3 7-16)/( 7 -4) 2 



D(q,z) 



N (7-4) 5 



q 2 2 7 - 9 6 ( 9 2 ) 
7~^T 



37- 16 



V e F q J 2(7 - 4) 2 



and the polarizability x of the ideal ID Fermi gas with 
renormalized mass is given by the relation 



,«») 



(«,*) = £ 



nk+q/2 - n k -q/2 

hz — h 2 kq/m* 



(29) 



In the thermodynamic limit wc find the real and imag- 
inary parts of x(°)(q,(jj + is) = x{ (q>w) + ixi (q^m) 
using the relation l/(x + is) = P(l/x) — iird(x) 



Xi 



*2 (<?. W ) 



TVm* 
2H 2 qk F 

Nm* 
2H 2 qk F 



dk 



nk+q_ — n k + q+ 



w(n q _ -n q+ ), 



(30) 
(31) 



where P means the Cauchy principal value and wc defined 



q± 



Hq 



4 



(32) 



At zero temperature the occupation numbers n k define 
the Fermi step function and we arrive at the simple ana- 
lytic expressions 



X ( i\q,uj) = 

X 2 0?> W ) = 



Nm* 
2H 2 qk F 

Nirm* 



In 



^-(q) — u 2 
J ±1, ui- ^5 ±cj ^ w+, 
2H 2 qk F I 0, otherwise. 



The dispersion relations 

w±(q) = H\2k F q±q 2 \/(2m*) 



(33) 
(34) 

(35) 



border the continuum part of the accessible excitation 
spectrum made up from HF quasiparticle-quasihole ex- 
citations i|16fl . a s shown in Fig. [21 The two branches 
u>± (q) thus approximate the two branches of elementary 
excitations introduced by Lieb Q as type I and type II 
excitations, respectively. 

In accordance with the exact results, both branches 
share the same slope at the origin and give rise to 
a single speed of sound at zero temperature given by 
vt = du±/dk = hk F (l — 47~ 1 )/m. This value is the cor- 
rect first order expansion [2j of vt for large 7, consistent 
with Eq. I|22(l . Note that the usual Bogoliubov pertur- 
bation theory |4l| for weakly interacting bosons gives a 
similar expansion of vt for small 7 and the type I ex- 
citation branch. Type II excitations are not described 
with Bogoliubov theory. The dispersion curves ui± (q) of 
Eq. 135|) differ from the free Fermi gas (TG gas) values 
only by the renormalization of the mass, which already 
takes place in the HF single-particle energies. 



B. Dynamic structure factor 

The DSF S(q, ui) is the Fourier transform of the 
density-density correlation function [3^, |3^, |4jj an d ex- 
presses the probability to excite a particular excited state 
through a density perturbation 

S(q, uj) = Z- 1 e-P E ™\(m\p q \n}\ 2 5(Kw - E n + E m ), 

(36) 

where p q = J2i ex P( — iqxi) is the Fourier component of 
the density operator, Z = J2 n exp(—(3E n ) is the parti- 
tion function. 

The DSF is related to the dynamic polarizability by 



X(q,uj + is) 



\W (37) 
to' 2 -(u + is) 2 K ' 



7 




FIG. 2: (Color online) The excitation spectrum and the DSF 
at 7 = 13 for various temperatures. The upper and lower 
thin (blue) lines show the dispersions u+(q) and of 
Eq. (13511 . respectively, limiting the elementary excitations of 
the Lieb-Liniger model at T = 0. The dimensionless value 
of the rescaled DSF S(q, w)gei?/(fcpAT) from Eqs. (EOl and 
I43H is shown in shades of grey between zero (white) and 1.0 
(black). Here kp = Tin and ep = h 2 kp/(2m). The dotted 
(red) line indicates a 5- function contribution at u>o(q). At 
non-zero temperatures the 8- function contribution washes out 
and becomes a part of the continuum. 



or, equivalcntly |40j| . by 



S(q,w) 



Imx(q,oj + is) 



7r[l — exp(— fltku)] ' 
which gives at zero temperature 



(38) 



S(q, 



lmx(q, oj + ie)/n, oj > 0, 
0, oj < 0. 



(39) 



1. Zero temperature 

The DSF at zero temperature can be obtained from 
Eqs. (J3SJ) and which result in 



S(q,u) 



7r(l - 4 7 -!) 



B 



+ 



K o) y 



+ 6[(J-w (q)}A(q)/h, 



(40) 



with X12 given by the zero temperature expressions i|33|) 
and (|34|l . A grey scale plot of this result is shown in 
Fig. H The DSF of Eq. gJJ has two contributions. 
The first part is continuous and takes nonzero (and pos- 
itive) values only for oj_ < oj < oj +1 which is also the 
region where particle-hole excitations on the HF level 
are present. The second part is a discrete branch with 
strength A{q) and located at oj = u>o{q), outside the re- 
gion of the discrete contribution. As we will discuss in 
detail below, the discrete part is exponentially suppressed 
for small 7 -1 and should be understood as an artefact of 
the RPA approximation. 

Due to a logarithmic singularity in Xi\ the DSF van- 
ishes on the dispersion curves oj±(q). For the TG gas at 
7 — > 00, the value of the DSF within these limits is inde- 
pendent of oj and takes the value of Nmj (2irh 2 qn). The 
energy-dependence in the RPA for finite 7~* is shown 
in Figs. [3 and |31 In particular, we see that the umk- 
lapp excitations at q = 2fcp and small oj, which prohibit 
superfluidity of the TG gas, are suppressed for finite 7. 
We find that S(2kp,oj) in the RPA approaches zero as 



1/ In (Hoj/ef), in contrast to the results of Refs. (42ll4 
which predict a power-law dependence on u for finite 7 
based on a pseudoparticle-opcrator approach. 

In the RPA the enhancement of Bogoliubov-like exci- 
tations is seen as a strong and narrow peak of the DSF 
in the RPA near oj+ at large momenta in Fig. |3J At 
finite gamma and for small momenta q < nn/2, how- 
ever, the RPA predicts a peak near w_, in contrast to 
the first-order result. Whether this effect is real or an 
artefact of the RPA is not obvious and may be decided 
by more accurate calculations or experiments. Spurious 
higher order terms in the RPA and an improved approx- 
imation scheme have been discussed in Ref. 01- O n the 
other hand, Roth and Burnett have recently observed a 
qualitatively similar effect in numerical calculations of 
the DSF of the Bose-Hubbard model 0. 

The RPA result may be expanded in I/7, which is con- 
sistent with direct perturbation theory up to first order. 
This yields for lj_ ^ oj $J oj+ 



S(q,oj) 



N 



I + 87- 
4g 



ln/(g,^) 
2 7 



+ 0{^ 2 ) (41) 



with f(q,oj) = \(oj 2 — oj 2 _)/{oj 2 _ — oj 2 )\. However, this 
first order expansion can assume negative values as seen 
in Fig. El although the DSF, given by Eq. (J3BJ), should 





FIG. 3: (Color online) The DSF S(q,u) ef/N at interaction strength 7 = 13.33 as a function of ui at q — 0.4fcp and q = 2/cf 
for various temperatures. The solid (black) line shows the RPA result 1401 or 143H . the dashed (red) line shows the first order 
expansion in 7 -1 141H or 14411 . and the dotted (blue) line shows the DSF of the TG limit (7 = 00) for comparison. The circular 
(violet) dots on the x-axis denote u)±. The first-order expansion always has a divergence to ±00 near u± at T = 0, respectively, 
and the maximum is shifted to ui+ at finite T. The RPA result, however, shows an enhancement of low-energy excitation near 
uj- at small q. Note the unphysical negative values of the DSF in the first order expansion near in particular for the 
umklapp excitations at q = 2&f and u> close to zero. 



be strictly non-negative, a property that our RPA result 
(|40|l fulfills. Close to the first order expansion has a 
logarithmic singularity tending to +00, which may be a 
precursor of the dominance of Bogoliubov-like excitations 
in the DSF at small 7. In the TG limit 7- 1 = 0, the DSF 
becomes discontinuous with respect to uj at u>± because 
the DSF of the TG gas is a step function [see Eq. . As 
a consequence, the first order approximation (|41|l cannot 
be good for arbitrary values of q and uj but diverges in 
vicinity of ui± due to the slow convergence of perturbation 
theory close to the point of discontinuity. There is no 
formal problem here since the expression (|41|) remains 
positive if for any given finite value of q and uj 7^ uj- a 
large enough value of gamma is chosen. 



Finally we discuss the (5-function part of the DSF (|4l)|l . 
This contribution relating to discrete excitations of col- 
lective character in the time-dependent HF scheme lies 
outside of the continuum part and comes from possi- 
ble zeros in the denominator of x(g, uj + is). It is de- 
termined by the solution Wrj(s) of the transcendental 
equation B = — D(q, uj)Xi (OS in conjunction with 

X2 few) = 0. We have solved this equation in vari- 
ous limits and found that at most one solution for u>a(q) 
exists. The strength A(q) is given by the residue of the 
polarizability at the pole z$ = u)o(q). After small algebra, 



we derive from Eq. (12811 



A(q) = N 



(7~4) 3 



4(37-4) huj {q) 1 + 16?7 2 V 



(42) 



where rj = q/kv and 



hn = 



7(7 - 6) 2 [ In |£ - (V - 2) 2 | - In |£ - („ + 2) 2 |] ~' 



3 7 -4 



[£-(r/-2) 2 ]K-(77 + 2) 2 



with £ = [fi Wo ( g V£F] V/N(7 ~ 4)] 2 . 

Numerical values for A{q)ujo{q) are shown at finite 7 
in Fig. 21 For small q we find a ^-function contribution 
at u>o{q) < uj- whereas for large q there is a discrete 
contribution at uo(q) > uj+ (see Fig.|2Ji). In the limit q — > 
00 at finite 7, the <5-part completely determines the DSF 
as the continuum part vanishes; asymptotically A ~ N , 
and wo — hq 2 / (2m) becomes the free particle dispersion, 
reminiscent of the DSF for the weakly interacting Bosc 
gas at large momentum in Bogoliubov theory |4l| . 

For small 7 -1 , the strength A(q) ~ 2Nj exp(—jq/kp) 
is exponentially suppressed and possible solutions are 
close to the dispersion branches uj± with \ujo — uj± \ oc 
cxp(— -fq/kp). Due to this proximity of the discrete and 
continuous parts and expected smearing of discrete con- 
tributions by interactions beyond the RPA, we may con- 
jecture that the 6-i unction should be seen as part of 
the continuum, enhancing contributions near the bor- 
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FIG. 4: (Color online) Contour plot of A(q) 2mu) (q) / (hNq 2 ), 
which shows the importance of the <5-function contribution. 
Values are given in ten contours between (white) and 1 
(black). Between the thick (red) lines there is no discrete 
contribution because \ given by Eq. 128H has no poles. Left 
of this region, there is a discrete part with uiq < w_ and right 
of it there is one with too > w+. Above the dashed (green) line 
the RPA breaks down due to an instability of the HF ground 
state as discussed in Sec. Mil 



frequency, the DSF decays exponentially in accordance 
with Eq. Similar to the case of zero temperature, 

the enhancement of excitations still takes place close to 
ui+ for q > 7m/ 2 and w_ for q < nn/2 at small values of 
temperature T < 0.5£f- 

To the first order in 7 -1 , we have 



'</+ 



1 — exp(— pfiw) 



I + 87- 



■Ip / dfc nk+ «+ - 
27 i k 



rik+ q _ 



(44) 



This linear approximation fails in vicinity of the umklapp 
excitation q — 2fcp and to = 0, yielding unphysically neg- 
ative values of the DSF contrary to the obtained RPA 
expression see Fig.|3] 



3. Sum rules for the DSF 



der. Moreover, at finite temperatures there is no 5- 
function contribution even within the RPA, as we discuss 
in Sec. IIVB2l below. Indeed, we know from the exact so- 
lutions |2j that the energy spectrum is continuous. 

The RPA polarizability l|28|) is a retarded Green's func- 
tion and thus has to be analytic in the upper half com- 
plex plane |38l.l39j]. At zero temperature, the analyticity 
breaks down above the dashed (green) line in Fig.0] The 
instability of the RPA results from the instability of the 
HF approximation |46[ and arises exactly at the critical 
value of 7 when the isothermal speed of sound equals to 
zero, see Fig. ^ 



2. Finite temperatures 

At finite temperatures we obtain the DSF by means of 
Eqs. (H| and ® 



S(q,uj) 



7T(1 -47- 1 ) 



(43) 

where the real and imaginary parts of the polarizability 
x[ 2 are given by Eqs. (|3Tj|) and (|3*T|) . respectively. The 
DSF is shown in Figs. |21 and [3] The main effect of fi- 
nite temperature is a smoothing of the zero-temperature 
features. The (5-function contribution to the DSF dis- 
appears, since x 2 7^ an d thus the denominator of 
Eq. I|43|) does not vanish for lu ^ 0. It is absorbed by 
the continuum part of the DSF. 

At finite temperatures, the non-vanishing contribu- 
tions of the DSF spread considerably beyond the particle- 
hole excitation spectrum limited by a->_ and uj + because 
the DSF no longer probes the ground state but a ther- 
mal ensemble [see Eq. I|36[) ]. For negative values of the 



Sum rules for the DSF are an important test for check- 
ing the validity of the obtained expressions. In particular, 
the /-sum rule [H, US E3| 

mi = h 2 J uS{q, tu)duj = Nh 2 q 2 /{2m) (45) 

should be fulfilled to all orders in 7 _1 within the 
RPA 0| . We have verified it by numerical integra- 
tion and found excellent agreement at finite values of 
7 for both the zero-temperature DSF l|40|) and the fi- 
nite temperature expression (|43p . The /-sum rule can 
also be verified analytically from the large u> asymptotics 
X(q,w) ~ -2m 1 /(Huj) 2 using Eqs. ancl (EH), assum- 
ing that x is analytic as a function of u> in the upper half 
complex plane. 

The sum rule for the isothermal compressibility |40l | 



limP 

9^0 



S(q,c 



-duj = 



N (dn\ 



2n \dfiJ t 



(46) 



holds also to all orders in 7 -1 , which can be checked 
analytically. Indeed, by comparing Eq. I|28|) with the HF 
isothermal compressibility discussed in Sec. lIIII wc derive 



I n ■ ^ N ( dn 

hm x{q, + is) = — — 
q->o n \dfi/T 



(47) 



Then Eq. 1|46|) is a direct consequence of the dispersion 
relation (|37() at ui = 0. 



4- Consequences for superfluidity 

As we have argued in Ref. [3l]], the value of the DSF 
near the umklapp excitations at uj = and q = 2kp is 
relevant for the phenomenon of superfluidity according 
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to the Landau criterion. A finite value of S(q,u>) will 
prohibit persistent currents, as spontaneous excitations 
initiated by infinitesimal perturbations would be able to 
dissipate the translational kinetic energy stored in the 
current. Our finite temperature results clearly show that 
S(q = 2kp,uj = 0) at a given value of T > will al- 
ways be positive and finite for large enough 7 in both 
the RPA expression l|43|) and the first order expansion 
We thus conclude that there is no superfluidity in 
the large-7 regime at finite temperatures, in accordance 
with Popov's analysis made many years ago. The 
question of efficient suppression of the DSF in the vicin- 
ity of the umklapp excitation at T = cannot be fully 
answered within the present approach due to the nonreg- 
ularity of S(q,u>) at this point, as discussed in Ref. [3ll |. 
and will be left to future investigations. 



C. Static structure factor and pair distribution 
function 

The static structure factor S(q) [38l Isfj l4(1 | is a func- 
tion of momentum only and is obtained by integrating 
the DSF over the frequency 



S(q) 



(PqP, 



N 



h 

N 



duj S(q, to) 



(48) 



The results of numerical integration of the DSF in the 
RPA are plotted in Fig. 

The static structure factor contains information 
about the static correlation properties of a system 
and directly relates to the pair distribution function 
or the normalized density-density correlator g(x) = 
(*+(x)* t (0)«'(0)4'(x))/n 2 by the equation 



1 



(49) 



At small momenta S{q) can be related to the isother- 
mal compressibility because the main contribution into 
the integral (|48J) comes from the "classical" region H(3oj <C 
1 H3- We derive from Eqs. ©, I03l, and gHJ) 



q->0 n V o\l' T TOWJ, 



(50) 



This relation of the structure factor S(q) to the speed 
of sound i>t at small momentum implies that S(q) in the 
RPA is exact up to first order in 7 _1 and is overestimated 
at finite 7 as seen from the results for vt in Fig. ^ 



1. Zero temperature 

We can obtain the static structure factor (|48|l from 
Eq. l(4*T|) to the first order 

S(q)=SW(q) +j- 1 S 1 (q) +0(^ 2 ). (51) 



Here S 1 ' ' denotes the static structure factor for the ideal 
ID Fermi gas 

./M/PM-IK*. (52) 

\ 1, |«| > 2k r , 

and the function Si (q) takes the form 

Si(q) = M Hv) -\V~2\ In \ V - 2| - \r, + 2| In \ v + 2|] 
+ 4S {0 \q) (53) 

with the dimcnsionless wave vector 77 = q/kp and the 
function 



r(r?) = 



41n2, \ V \^2, 
2|r ? |ln| 7? |, \ V \ 2 2. 



(54) 



The obtained correction S±(q) is continuous and has 
the asymptotics Si(q) ~ — 8/(3r/ 2 ) + 0{l/rf) when ?] — > 
00 and Si(q) ~ 2|r?| - \rj\ 3 /2 + 0(r? 5 ) when rj -> 0. The 
latter asymptotics gives us a possibility to determine the 
coupling constant 7 experimentally from the phonon part 
of the static structure factor (|51|l for q < fcp 



S(q) 



2fo 



4 

1 



(55) 



The static structure factor at small q is related to the 
sound velocity by S(q) ~ H\q\/(2mvT)- It is easily 
seen that our result (|55|l is consistent with the sound 
velocity of Eq. (|22)> . 

Figure |5l shows the static structure factor in the full 
RPA and its first order expansion (|51|) at 7 = 13.3 and 
the TG limit for comparison. The first order re- 
sult shows a cusp which is an artefact of the first order 
expansion. 

Using Eq. (|49|l in conjunction with relations (|51|) - (|54|l . 

we can represent our result for the pair distribution func- 
tion in the form 



g{x)= 1- 



. 2 

sin z 



2n d sin 2 z 



4 sin 2 z 



2d_ 

7 dz 



sin z 



7 dz z 2 
dr] sin(?7z) In 



7 z 

1+7? 



1 



(56) 



where z = kpx = mix. It follows from this equation 
that g(x — 0) vanishes not only in the TG limit but 
also in first order of 7 , which is consistent with the 
results of Refs. and the HF expression H60|) below, 

indicating once more the validity of our results. The 
physically correct limit g(x) — » 1 for x — > 00 is fulfilled 
due to Eq . 1491) . A similar expression for g(x) was derived 
in Ref. [15| for the large distance asymptotics. To our 
knowledge, Eq. (|56|) shows for the first time the full x 
dependence of g{x) up to first order in 7 _1 . 



1 1 




q/k F q/k F q/kp 



FIG. 5: (Color online) The static structure factor S(q) as a function of momentum for 7 = 13.3 and various temperatures. 
The solid (black) line shows the RPA result obtained by numerical integration from Eqs. 14U1 or 14311 . The dashed (red) line 
shows the first order expansion 1511 in j~ 1 with Eqs. I52H and 1531 at zero temperature or Eqs. 1571 and 1581 1 at non-zero 
temperatures. The dotted (blue) line shows the static structure factor in the TG limit (7 = 00) 1521 or 157H for comparison. 
Note the unphysical cusp at q = 2kp in the first-order expansion of the static structure factor at zero temperature and some 
traces of it at small temperature. 



2. Finite temperatures 

The first-order approximation for the static struc- 
ture factor at finite temperatures is obtained by using 
Eqs. O and 0- The result takes the form of Eq. JST} 
but with the TG, or free Fermi, static structure factor 



2hkpq 

and with the function S±(q) 



1 — exp(— (3fruj) 



(57) 



Sl{q)= lk 



dcu 



'<i - 



1 — exp(— (3Huj) 



P / dfc 



■4S<°>(9), 



(58) 



where q± is given by Eq. (|32f> at 7 _1 = 0. The finite 
temperature results are plotted in Fig. El One can see 
an unphysical behaviour of the first order approximation 
near q = 2k-p, in contrast to the full RPA result. The 
small momentum limits for S(q) arc determined by the 
speed of sound «t through Eq. J5UJ. In the plot on the 
right hand side, the RPA overestimates S(q) at small q 
as a result of the deviations of the speed of sound as seen 
in Fig.Hl 



3. Limits of validity 

When the interaction is proportional to a small param- 
eter, the RPA method is applicable and yields correct val- 
ues of the DSF at least up to the first order in this param- 
eter [32l . This implies the validity of the obtained 
expressions for the polarizability, the dynamic and static 
structure factors, and the pair distribution function up to 
the first order in 7 . Smallness of the inverse coupling 
constant means, in particular, small values of the pair 
distribution function in the contact point: g(x = 0) < 1, 
which is the TG regime by definition. 

A classification of different regimes in the ID Bose gas 
for arbitrary temperatures was given in Rcf. |20|, As it 



was mentioned above, it is possible to obtain the values of 
g(x) at x = from the exact solution of the Lieb-Linigcr 
model with the help of the Hellmann-Feynman theorem. 
The TG regime is realized [2(| when 



7 > max(l, y/T/e F ), 



(59) 



which gives also the criterion of validity of the RPA re- 
sults. We can derive this criterion within the HF ap- 
proach of Sec. IIIII Indeed, we obtain the following re- 
sult for the pair distribution function by applying the 
Hellmann-Feynman theorem to the HF grand potential 



g(x = 0) 



Air 2 (fc 2 ) 



(60) 



Using the low-temperature expansion of the average mo- 
mentum (fc 2 ) = (fc|/3)[l + 7r 2 T 2 /(4e 2 ) + ■ ■ •] and 
the high-temperature expansion (fc 2 ) = fc|T/[2ep(l — 
47 _1 )] + • • • , we arrive at the above mentioned restriction 
on 7. 

The validity of the RPA requires, in particular, the sta- 
bility of the HF solutions [4g. Thus our results are ap- 
plicable in practice for 7 > 10, see discussion in Sec. 11111 

Note that the HF expression l|60(l yields the correct 
value of the pair distribution function only at x = but 
up to the second order in 7 _1 . This is due to validity of 
the HF approximation in the first order in 1/7; hence, the 
derivative with respect to 7 of the first-order correction 
for the grand potential gives the correct value of g(x = 0), 
proportional to I/7 2 . By contrast, the RPA expression 
(|56[) and its finite temperature generalization yield the 
values of g{x) for arbitrary x but guarantee validity only 
up to the first order. For this reason, the numerical values 
of g(x = 0) obtained with Eqs. (15jl. (Igjl. and gHJ differ 
from those of Eq. IG(JI) . 
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V. CONCLUSION 

We have derived variational approximations for the dy- 
namic polarizability and related two-particle correlation 
functions of the one-dimensional Bosc gas, extending our 
previous results |3l| to finite temperatures. The approx- 
imations are good for strong interactions and yield ex- 
pansions valid to first oder in 1/7, which had not been 
available previously. We have carefully checked the con- 
sistency with known limits and sum rules and analyzed 
the limits of validity of the derived equations. Due to 
the Bose-Fermi duality, our results are equally applica- 
ble for strongly interacting bosons as well as for weakly- 
interacting spinlcss fermions. Our result for the DSF in- 
dicates a dramatic departure from the TG limit already 
for very small values of I/7 by enhancing Bogoliubov- 
like excitations and by suppressing umklapp excitation, 
which are the main obstacle to observing superfluid-likc 
response in the ID Bose gas. However, we find that su- 



perfluidity at finite temperatures is strictly prohibited in 
the large-7 regime as umklapp excitations are always as- 
sociated with a finite probability. Finite temperature ef- 
fects generally are found to smear out the sharp features 
of the zero temperature correlation functions. Neverthe- 
less, at a level of 10% of the Fermi temperature £f, the 
main effects should be well observable in experiments. 

Our results also establish the usefulness and validity 
of the fermionic pseudopotcntials (jHJ) and |(2J and the 
variational Hartrce-Fock approximation and RPA. The 
method can easily be extended to further studies in the 
large-7 regime by including the effects of harmonic or 
periodic external potentials or by studying nonlinear re- 
sponse properties. Furthermore, the acquired knowledge 
of the dynamic density correlations will be useful for con- 
structing an accurate time-dependent density functional 
theory, extending the approach of Ref. • 

The authors are grateful to Sungyun Kim and Rashid 
Nazmitdinov for useful remarks. 
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